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ABSTRACT 

The acquisition of late-time imaging is an important step in the analysis of pre-explosion 
observations of the progenitors of supernovae. We present late-time HST ACS WFC ob- 
servations of the sites of five Type IIP SNe: 1999ev, 2003gd, 2004A, 2005cs and 2006my 
Observations were conducted using the F435W, F555W and F814W filters. We confirm the 
progenitor identifications for SNe 2003gd, 2004A and 2005cs, through their disappearance. 
We find that a source previously excluded as being the progenitor of SN 2006my has now 
disappeared. The late-time observations of the site of SN 1999ev cast significant doubt over 
the nature of the source previously identified as the progenitor in pre-explosion WFPC2 im- 
ages. The use of image subtraction techniques yields improved precision over photometry 
conducted on just the pre-explosion images alone. In particular, we note the increased depth 
of detection limits derived on pre-explosion frames in conjunction with late-time images. We 
use SED fitting techniques to explore the efifect of different reddening components towards the 
progenitors. For SNe 2003gd and 2005cs, the pre-explosion observations are sufficiently con- 
straining that only limited amounts of dust (either interstellar or circumstellar) are permitted. 
Assuming only a Galactic reddening law, we determine the initial masses for the progenitors 
of SNe 2003gd, 2004A, 2005cs and 2006my of 8.4 ± 2.0, 12.0 ± 2.1, 9.5!|^ and 9.8 ± 1.1 Mq, 
respectively. 

Key words: stars : evolution - supernovae : general - supemovae : individual : 1999ev - 
supernovae : individual : 2003gd - supernovae : individual : 2004A - supemovae : individual 
: 2005cs - supernovae : individual : 2006my 



1 INTRODUCTION 

AH stars with initial masses > 8Mo are expected to end their lives 
as core-collapse supemovae (CCSNe). In the last decade, the di- 
rect observation of the progenitors of CCSNe in fortuitous pre- 
explosion imaging, has become an integral step in the study of all 
nearby events (for a review see Smartt 2009). 

The majority of the success in the actual detection of progeni- 
tors has been for the Red Supergiant (RSG) precursors to hydrogen- 
rich Type II Plateau (IIP) SNe. The ensemble of progenitor de- 
tections and detection limits for these RSGs, however, presented a 
conflict with the theoretical expectation from stellar evolution mod- 
els (Smartt et al. 2009). The "Red Supergiant Problem" refers to 
the apparent absence of progenitors of Type IIP SNe with mass 
> 17Mo, whereas the theoretical expectation is that stars with 
masses up to 25 - 35Mo should end their lives as RSGs. Recently 
Smith et al. (2011) and Walmswell & Eldridge (2012) considered 
the role of circumstellar dust, which is otherwise not probed by the 
resulting SN or the surrounding stellar population, as a possible so- 
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lution to the RSG problem. In particular, Walmswell & Eldridge 

established that the amount of circumstellar dust, and hence red- 
dening, is larger for higher luminosity RSGs, that arise at the upper 
end of the mass range for stars to explode as Type IIP SNe. In the 
case of SN 2012aw, Eraser et al. (2012) and Van Dyk et al. (2012) 
observed a significant difference between the reddening determined 
towards the progenitor and the reddening inferred towards the sub- 
sequent SN; suggesting a significant amount of dust was associated 
with the progenitor and that it was destroyed in the SN explosion, 
such that it could not be measured post-explosion. Kochanek et al. 
(2012) and Van Dyk et al. (2012) also considered the additional 
possibility that the nature of the reddening, due to a circumstellar 
dust component, was different to ordinary reddening laws appro- 
priate for the interstellar medium (e.g. Cardelli et al. 1989). 

A further issue with the RSG problem is that the maximum 
mass is poorly constrained, due to the paucity of progenitors with 
high inferred masses. This may reflect either a real deflcit of high 
mass progenitors (i.e. that the upper mass limit for stars to explode 
as Type IIP SNe is low) or may be due to low numbers of high mass 
progenitors according to the Initial Mass Function (IME). The exact 
articulation and quantification of the RSG problem is compounded 
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by the reliance on, generally, poor fortuitous pre-explosion images, 
limited detections in pre-explosion images in multiple filters, and 
ambiguity as to whether the object observed at the SN position is 
indeed the progenitor object. 

Maund & Smartt (2009) demonstrated that it is possible, us- 
ing deep late-time images that are tailored (in a way that the pre- 
explosion images, by their fortuitous nature, cannot be) to observe 
the site of a given SN and overcome the severe limitations of anal- 
ysis on just the pre-explosion images alone. Using image subtrac- 
tion techniques, such as isis (Alard & Lupton 1998; Alard 2000), 
the late-time images can be used as templates, to both demonstrate 
the disappearance of the progenitor candidate (confirm the original 
identification) and conduct precise photometry (without contami- 
nation from nearby and underlying objects). 

Here we present a new analysis of five previously identi- 
fied progenitor candidates using newly acquired late-time HST 
Advance Camera for Surveys (ACS) Wide Field Channel (WFC) 
imaging. These progenitors are for SNe 1999ev, 2003gd, 2004A, 
2005cs and 2006my, and details of these SNe are presented in Ta- 
bid. 



2 OBSERVATIONS 

The analysis of the pre-explosion HST observations of the sites 
of the five target SN progenitors has been previously presented by 
Maund & Smartt (2005) (1999ev); Van Dyk et al. (2003b), Smartt 
et al. (2004) and (Maund & Smartt 2009) (2003gd); Hendry et al. 
(2006) (2004A); Maund et al. (2005) and Li et al. (2006a) (2005cs); 
and Li et al. (2007), Leonard et al. (2008) and Crockett et al. (2011) 
(2006my). The pre- and post-explosion and late-time observations 
of the sites of the five target Type IIP SN progenitors are presented 
in Table 2. 



2.1 Pre-explosion observations 

Due to the availability of more recent and appropriate calibrations, 
new versions of the pre-explosion WFPC2 and ACS data were re- 
trieved from the HST archive' . The pre-explosion observations ac- 
quired with WFPC2 were drizzled together following the standard 
procedure^. These images were produced for the purposes of image 
subtraction (see Section 2.6). In parallel, a separate "reduction" and 
photometric analysis was conducted on the same data using HST- 
phot (Dolphin 2000b). For SN 2005cs, pre-explosion images were 
acquired with the Advanced Camera for Surveys (ACS) Wide-Field 
Channel (WFC). These images were acquired using a four-point 
box dither pattern. Alignment between the images was checked us- 
ing the PyRAF task tweakshifts, and the images were combined 
using multidrizzle. Due to the half-integer pixel shifts of the dither 
pattern it was possible to enhance the spatial sampling of the final 
image, providing a final pixel scale of 0.035". This is larger than 
exact half-sampling (0.025" px '), but is used to match the sam- 
pling achieved for the late-time ACS WFC images of SN2005cs 
(see Section 2.3). 



http://archive.stsci.edu 
^ http://www.stsci.edu/hst/wfpc2/analysis/WFPC2_drizzle.html 



2.2 Post-explosion observations 

For the purposes of this study, the principal interest in the immedi- 
ate post-explosion images (acquired up to 3 years post-explosion) 
was to provide a position for the SN relative to the surrounding 
stars, such that the SN position could be identified on the pre- 
explosion and late-time frames through differential astrometry. The 
"reduction" procedure for these images was the same as outlined 
for the pre-explosion observations (see Section 2.1). 

2.3 Late-time ACS observations 

Late-time observations of the sites of four of the target Type IIP 
SNe were acquired using the HST ACS WFC 1 for the program 
GO-11675 (PI Maund). The WFC chip was windowed to an array 
of Ik X Ik pixels to reduce the readout time and mitigate the role of 
Charge Transfer Inefiiciency (CTI). Observations were conducted 
in three filters FA35W, F555W and FilAW. Importantly, for two 
SNe (2003gd and 2004A) the use of the F555W^ filter in these 
late-time observations does not match the pre-explosion images ac- 
quired with the wider F606 W filter. In the pre-explosion frames, the 
majority of the flux from the progenitor at these wavelengths is rep- 
resentative of the continuum (requiring a small F5551V - F6061V 
colour correction). At late-times, however, the F606iy filter en- 
compasses the wavelength of Ha, which is a characteristic emis- 
sion feature of late-time Type IIP SN spectra. A more appropriate 
comparison between before and after continuum fluxes is, there- 
fore, achieved with the F555W filter although, as discussed below, 
corrections for the slightly different filter ttansmission functions 
also need to be considered. Each of these late-time observations is 
composed of four separate sub-exposures acquired in a 4-point box 
dither pattern. This arrangement was used to permit the acquisition 
of better spatial sampling of the point-spread function (PSF) and 
for removal of fixed hot-pixel features. The PyRAF task multidriz- 
zle was used to drizzle each of the sub-exposures (for a given filter) 
together. The task tweakshifts was used to fine-tune the aUgnment 
between each of the sub-exposures prior to drizzling. Although the 
4-point box dither pattern employs half-integer pixel shifts, poten- 
tially giving an improvement of spatial sampling by a factor of 2, 
the presence of aUasing (alternating bands across the frames) pro- 
hibited reaching this final pixel-scale. This phenomenon was rela- 
tively insensitive to the choice of the drizzling kernel. Instead the 
final pixel scale is greater than one half of the original 0.05" pixel 
scales of ACS/WFC (0.035" px~'). 

For SN 2006my, observed for program GO-12282 (RI. D. 
Leonard), the late-time images were only acquired in two bands 
(F555W and F814W). For each filter, two exposures were ac- 
quired, for the rejection of cosmic rays, but at the same pointing 
such that no spatial resampling was possible. These images for each 
filter were combined using multidrizzle, but with the output images 
having ttie origmal ACS/WFC pixel scale (0.05" px"'). 

2.4 Geometric Transformations 

A series of transformations were calculated for the complete 
datasets to determine the positions for objects in a common ref- 
erence frame. Geometric transformations were calculated between 
the pre-, post-explosion and late-time F555W images (or, if un- 
available, F814iy images) using the iraf task geomap, assuming 
only simple offsets, rotations and scalings. For the data at a given 
epoch, shifts between images with other filters and the correspond- 
ing reference F555W image were calculated by cross-correlating 
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Table 1. Details of the five target Type II SNe, with progenitor candidates identified in pre-explosion HST images. 



Supernova 


Type 


Host Galaxy 


Distance' 




[0/Hf 








(Mpc) 


(dex) 


1999ev 


II(P) 


NGC 4274 


15.1 ±2.6 


0.020 


8.5 


2003gd 


IIP 


NGC 628 (M74) 


9.3 ± 1.8 


0.061 


8.4 


2004A 


IIP 


NGC 6207 


20.3 ± 3.4 


0.014 


8.3 


2005CS 


IIP 


NGC 5194 (M51 


8.4 ± 1.0 


0.031 


8.7 


2006my 


IIP 


NGC 4651 


22.3 ± 2.6 


0.024 


8.7 



' After Smartt et al. (2009). 

^ After Schlafly & Finkbeiner (201 1), as quoted by NED. 
^Smartt et al. (2009). 



these images with the F555W image as the reference frame. The 

cross-correlation was facilitated using the PyRAF task crosscor, 
with the corresponding shifts calculated using shiftfind. 

2.5 Photometry 

For WFPC2 pre-explosion images the HSTphot package (Dolphin 
2000b) was utilised to conduct PSF-fitting photometry of the input 
images. HSTphot provides the appropriate corrections for aperture 
size and charge transfer inefficiency, as well as tools for conducting 
artificial star tests. We note that HSTphot only provides aperture 
corrections to a final aperture size of 0.5". Using the corrections of 
Holtzman et al. (1995), we apply a term to correct to the photometry 
to an infinite aperture. 

For data acquired for program GO- 1 1675 (PI Maund), princi- 
pal photometry was conducted using iraf DAOphot on the final out- 
put drizzled images with the enhanced spatial sampUng. We used 
the latest zeropoints appropriate for ACS WFC^. Aperture correc- 
tions were calculated for each image to an aperture of 0.5", with a 
further correction to infinity adopted from Sirianni et al. (2005). 

A key concern for the fidelity of the derived photometry is the 
inefficiency of charge transfer (CTI) for charged coupled detectors 
on HST. As photometry was conducted on the drizzled, subsampled 
images, evaluating the magnitude of the CTI on the final photom- 
etry is non-trivial and is based on the position of a given star on 
the original undrizzled, distorted FLT images. Geometric transfor- 
mations were calculated between the final drizzled images and the 
FLT images, using 3rd order polynomials in x and y, using ge- 
omap. This approach was used, over using a simple pre-computed 
distortion table, as non-negligible shifts were found between the 
expected pointings in the dither pattern. The positions of stars on 
the output drizzled images were transformed to the corresponding 
locations on each of the four constituent FLT images. Following 
Annibali et al. (2008), we conducted small aperture (3px) photom- 
etry on the individual FLT images and used the measured flux and 
sky background values to calculate the magnitude loss due to CTI 
following the analytic prescription of Chiaberge et al. (2009). Prior 
to conducting aperture photometry, the FLT images were scaled 
with the corresponding Pixel Area Map For a given star it might 
be only possible to calculate the magnitude loss for 1 or 2 of the in- 
put FLT images, because of the relative positions of bad pixels or 
cosmic rays. An average CTI loss (as a magnitude) was determined 
over the four constituent FLJ images and appUed to the photometry 
derived from the drizzled images. 



In addition, photometry of the ACS images was also con- 
ducted using the DOLPHOT package ^. We utilised two implemen- 
tations of DOLPHOT for photometry of the ACS data: DOLphot 
with the ACS module on the distorted CRJ and FLJ frames (which 
we refer to as DOLPHOT/ACS) and DOLPHOT as a generic pho- 
tometry package for the distortion-corrected drizzled frames. We 
find excellent agreement between our DAOphot photometry and 
the photometry derived using HSTphot and DOLPHOT, within the 
limits of the photometric uncertainties. Similarly to HSTphot, the 
ACS photometry was corrected from a 0.5" to an infinite aperture, 
using the corrections tabulated by Sirianni et al. (2005) 

The data for SN 2006my were analysed separately using only 
the HSTphot and DOLPHOT packages for the WFPC2 and ACS 
data, respectively. 



2.6 Image Subtraction 

Image subtraction techniques were used to conduct template sub- 
traction of the late-time images from the pre-explosion images to: 
1) confirm the identities of the progenitors through disappearance; 
and 2) conduct optimal differential photometry, independently of 
the background, of the now absent progenitors. We utilised the 
ISIS v2.2 image subtraction package (Alard & Lupton 1998; Alard 
2000), which matches the PSFs of the input and template/reference 
images (as well as refining the alignment between the two images 
and scaling the flux levels). In addition, isis also provides automatic 
object detection and photometry on the difference images. We also 
tested our image subtractions using the HOTPANTS image subtrac- 
tion package*, and examined the resulting difference images using 
DAOphot. We found no systematic difference between the photom- 
etry of difference images calculated using the two packages, and 
for this study use photometry derived from difference images con- 
structed using ISIS 

For each SN the late-time images were, generally, used as the 
reference images. The late-time images were transformed to match 
the pre-explosion images using iraf task geotran. This ensured that 
we avoided resampling the lower quality pre-explosion images to 
match the superior late-time images. 

The evaluation of the systematic uncertainties was conducted 
by varying the key parameters in ISIS that principally affected the 
output photometry: the number and size of the stamps used for cal- 
culating the kernel, the degree of the kernel variation across the 
field and the degree of the background fitting function. A major 



^ http://www.stsci.edU/hst/acs/analysis/zeropoints/#tablestart 
http://www.stsci.edu/hst/acs/analysis/PAMS 



^ http://americano.dolphinsim.com/dolphot/ 

* http://www.astro.washington.edu/users/becker/hotpants.htmI 
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Table 2. HST observations of the sites of the five target Type II SNe. 







Date 


Instrument 


Filter 


Ljj\ IJIJ 3 Lli 

Time (s) 


Pi n n 1 Pi Y p 1 

Size (") 


PrrxTTflTTi 

JT 1 yJ^L (Xlll 


Pre-explosion 


U2JF0101/02T/03T 


SN 1999ev 
1995 Feb 5 WFPC2AVF2 


F555W 


280 


0.1 


5741' 


Post-explosion 


J8DT03010 
J8DT03020 
J8DT03030 


2001 Dec 31 
2001 Dec 31 
2001 Dec 31 


ACS/WFCl 
ACS/WFCl 
ACS/WFCl 


F555W 
F814W 
F435W 


450 
450 
400 


0.05 
0.05 
0.05 


9353^ 
9353 
9353 


Late-time 


JB4T01010 
JB4T01020 


2010 Nov 14 
2010 Nov 14 
2010 Nov 14 


ACS/WFCl 

ACS/WFCl 
ACS/WFCl 


F555W 

F814W 
F435W 


1368 

1408 
1608 


0.035 

0.035 
0.035 


11675' 
11675 
11675 


Pre-explosion 


U8IXCA01M/02M 
U8IXCY01M/02M/03M 


SN 2003gd 

2002 Aug 25 WFPC2AVF2 
2002 Aug 28 WFPC2AVF2 


F606W 
F606W 


1000 
2100 


0.1 
0.1 


9676^ 
9676 


Post-explosion 


J8NV01020 
J8NV01040 
J8NV01050 


2003 Aug 1 
2003 Aug 1 
2003 Aug 1 


ACS/HRC 
ACS/HRC 
ACS/HRC 


F435W 
F555W 
F814W 


2200 
1000 
1350 


0.025 
0.025 
0.025 


9733' 
9733 
9733 


Late-time 


JB4T02010 
JB4T02020 
TR4T02030 


2010 Nov 14 
2010 Nov 14 
2010 Nov 14 


ACS/WFC 
ACS/WFC 
ACS/WFC 


F555W 
F814W 
F435W 


1364 
1398 
1600 


0.035 
0.035 
0.035 


11675' 
11675 
11675 


Pre-explosion 


U6EAD001R/02R 
U6EAD003R/04R 


SN 2004A 

2001 Jul 2 WFPC2AVF3 
2001 Jul 2 WFPC2AVF3 


F814W 
F606W 


460 
460 


0.1 
0.1 


9042^* 
9042 


Post-explosion 


J8NV03010 
J8NV03020 
J8NV03030 


2004 Sep 23 
2004 Sep 23 
2004 Sep 23 


ACS/WFCl 
ACS/WFCl 
ACS/WFCl 


F435W 
F555W 
F814W 


1400 
1509 
1360 


0.05 
0.05 
0.05 


9733' 
9733 
9733 


Late-time 


JB4T03010 
JB4T03020 
TR4T03030 


2010 Sep 09 
2010 Sep 09 
2010 Sep 09 


ACS/WFC 
ACS/WFC 
ACS/WFC 


F555W 
F814W 
F435W 


1400 
1434 
1636 


0.035 
0.035 
0.035 


11675' 
11675 
11675 


Pre-explosion 


J97C5 
J97C5 
J97C5 
J97C5 


SN 2005CS 

2005 Jan 20-21 ACS/WFC 
2005 Jan 20-21 ACS/WFC 
2005 Jan 20-21 ACS/WFC 
2005 Jan 20-21 ACS/WFC 


F435W 
F555W 

F658N 
F814W 


2720 
1360 

2720 
1360 


0.035 
0.035 

0.035 
0.035 


10452"^ 
10452 
10452 
10452 


Post-explosion 


J9AR01011-31 


2005 Jul 24 


ACS/HRC 


F555W 


1944 


0.025 


11675 


Late-time 


JB4T04010 
JB4T04020 
JB4T04030 


2010 Jul 30 
2010 Jul 30 
2010 Jul 30 


ACS/WFCl 
ACS/WFCl 
ACS/WFCl 


F555W 
F814W 
F435W 


1460 
1494 
1696 


0.035 
0.035 
0.035 


11675 
11675 
11675 


Pre-explosion 


U2DT0901T/02T/03T 
U2DT0904T/05T/06T 


SN 2006my 

1994 May 20 WFPC2/WF2 
1994 May 20 WFPC2/WF2 


F555W 
F814W 


660 
660 


0.1 
0.1 


53751" 
5375 


Post-explosion 


U9OX030 1 M/02M/03M/04M 

U9OX0305M/06M 

U9OX0307M/08M 


26 Apr 2007 
26 Apr 2007 
26 Apr 2007 


WFPC2/PC 
WFPC2/PC 
WFPC2/PC 


F555W 
F814W 
F450W 


1200 
1200 
1400 


0.05 
0.05 
0.05 


10803" 
10803 
10803 


Late-time 


JBKSOlOlO 
JBKS01020 


21 Nov 2010 
21 Nov 2010 


ACS/WFC 
ACS/WFC 


F555W 
F814W 


1090 
1090 


0.05 
0.05 


12282'^ 
12282 




1 RI. J. Westphal 
5 RI. S.J. Smartt 
' P.I. M. Meixner 


2 RI. S.J. Smartt 
* RI. J. Rhoads 
RI. V. Rubin 


3 RI. J.R. Maund 
RI. S. Beckwith 
" RI. S.J. Smartt 


RI. S.J. Smartt 
* RI. FUippenko 
RI. D. Leonard 







Late- time view of progenitors of Type IIP supemovae 



concern for conducting image subtraction analysis using HST im- 
ages is tlie effect of tfie degree by wfiicfi tfie PSF is subsampled. 
Subsampling of the PSF means that the majority of the flux in 
WFPC2 images will fall in a single pixel, which may induce sys- 
tematic errors in the construction of the convolution kernel. In order 
to assess the systematic uncertainty associated with the degree of 
subsampling, iterations of the image subtraction routine were con- 
ducted with dififerent degrees of Gaussian smoothing applied to the 
input and/or reference images. We conservatively estimate that the 
total systematic uncertainties associated with the image subtraction 
process are ~ 5 - 10% of the flux observed in the difll'erence image 
(as discussed in Section 3; with the larger systematic uncertainty 
associated with the fainter progenitors, in addition to the commen- 
surate increase in the relative Poisson noise). The systematic uncer- 
tainties dominate over the Poisson noise, and we assume dominate 
over other noise sources, such as read noise, that are propagated to 
the difference images. 

As we are concerned with pre-explosion and late-time im- 
ages acquired with HST, the image subtraction process will involve 
CTI in both the input and reference images. There are two keys is- 
sues with the evaluating the effect of CTI on photometry that has 
been derived using difference imaging: 1) the CTI that affects the 
progenitor flux is in the flux system of the pre-explosion image, 
whereas the photometry of the progenitor from the difference im- 
age is found in the flux system of the late-time reference image; and 
2) the evaluation of CTI for a source using its photometry from the 
pre-explosion image explicitly undermines any increase in the pre- 
cision of photometry that might be derived using image subtraction 
techniques. 

For both WFPC2 and ACS/WFC the CTI is principally depen- 
dent on the flux of the object, the level of the nearby background, 
the position of the object on the chip (the number of charge trans- 
fers to be made to readout the electrons) and the date at which the 
observations were made. Rather than attempting to explicitly calcu- 
late the counts associated with the progenitor on the pre-explosion 
image, we instead use photometry of nearby or artificial stars, in 
the vicinity of the progenitor, to serve as proxies for the calculation 
of the CTI for the progenitor. As these stars are at approximately 
the same position on the chip, on similar backgrounds as the pro- 
genitor and observed at the same epoch, we can reduce the problem 
of determining the CTI to just the dependence on the brightness (or 
magnitude) of the progenitor The nearby stars (real or fake) sam- 
ple a range of brightnesses and simple expressions can be derived 
relating the magnitude of an object in the difference image to the 
CTT (in magnitudes) directly. 

We can consider the flux measured on the images, uncorrected 
for CTI, as /'. The flux corrected for CTI is then simply: 

f^flQ-OACTI 

Under the assumption that the CTI is a relatively small effect, we 
use the initial approximation that m « m' to derive the CTI in the 
pre-explosion frame using photometry, derived from difference im- 
ages, in the photometric system of the pre-explosion images. 

On the ACS images, we model the CTI as being effectively de- 
pendent on only two parameters: the brightness of the object and its 
y-position on the FLT images. We established, for the images con- 
sidered here, that the effect of nearby background inhomogeneities, 
around the progenitor positions, are well within the stated uncer- 
tainties of the CTI expressions. We consider the CTI for a given 
range of pixels to be approximately given by a power law depen- 
dent only on the magnitude of the object: 
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CTl(m) = 13m" (2) 

We evaluated the coefficients of Equation 2 for each pre-explosion 
ACS/WFC image (as the coefficients are dependent on the date on 
wliich the observations were made and the specific background at 
the progenitor position), using real and artificial stars for which 
the CTI had been evaluated using the equations of Chiaberge et al. 
(2009). Due to significant uncertainty in the expression used to de- 
termine the CTI, and its slowing varying nature with pixel position, 
it is possible to consider the CTI to be approximately fixed over 
~ 50 pixel ranges in y. 

We make a similar approximation for the WFPC2 observa- 
tions, using the CTI formulation presented by Dolphin (2000a)', 
derived using artificial stars generated using HSTphot. We find the 
average dependence of the CTI (over a 50 pixel range iny-position) 
to correspond to a second order polynomial: 

CTIim) = a+/3m + ym^ (3) 

For both WFPC2 and ACS observations, the dependence on 
m is relatively weak; the difference in its evaluation using m or 
m' is negligible. The importance of this approach is that it avoids 
specifically determining fluxes and sky background values from the 
pre-explosion images (which defeats both the purpose and preci- 
sion afforded by using image subtraction techniques to derive the 
photometry of the pre-explosion source). The zeropoint in the pho- 
tometric scale of the reference image Zg can be derived from pho- 
tometry rriR of reference stars, identified by isis. Using a package 
such as ISIS, of the reference objects in the reference image can 
be measured directly using aperture photometry and can be com- 
pared directly with the photometry of the same stars derived us- 
ing DAOphot, DOLPHOT or HSTphot; such fliat Zr contains not 
only the absolute zeropoint, but also all the relevant aperture cor- 
rections. The magnitude inj of the progenitor candidate can then be 
found directly. As the image subtraction procedure determines the 
difference in observed fluxes /', these fluxes must be further cor- 
rected for CTI derived on the input image (containing the progen- 
itor) using the scheme outlined above. Given the difference mea- 
sured from reference and input images, the final magnitude of the 
pre-explosion source is given as: 

= -2.5 log.o (/j) + Zr + CTI, (4) 

An additional source of systematic uncertainty is the differ- 
ences in the filter transmission functions between images used 
for image subtraction analysis. We note that, although some fil- 
ters are nominally identical, there may also be differences between 
the same filters used on different instruments. We used synthetic 
photometry of ATLAS9 (Castelli & Kurucz 2004) and MARCS 
(Gustafsson et al. 2008) model SEDs, using the total transmission 
(filter and instrument) functions, to determine the relative colour 
differences between the filter sets used here. For most combinations 
of filters, in particular between nominally identical HST filters, the 
colour difference as a function of temperature is < 0. 1 mags for low 
reddening (see Fig. 1). 

2.7 Non-detections and detection limits 

In previous studies (e.g. Maund & Smartt 2005; Crockett et al. 
2011), the derivation of the detection thresholds has been con- 
ducted using analytical expressions for the background and source 

' with updates from http : //purcell.as.arizona.edu/wfpc2_calib/ 
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Figure 1. Colour terms, between filter sets, as a function of temperature 
for ATLAS9 (heavy lines; Castelli & Kurucz 2004) and MARCS (grey line; 
Gustafsson et al. 2008) model SEDs (appropriate for supergiants with £(8- 
V) = 0; see Section 4). 

noise for an "ideal" observation. This approach, however, does not 
accurately reflect the way in which stars are actually detected in the 
photometry process, using iraf tasks such as DAOfind, and the ef- 
fect of crowding. We consider the insertion and attempted recovery 
of artificial stars to derive the detection threshold. 

Artificial stars were generated using the PSFs derived from the 
data themselves, with a randomly selected magnitude from a uni- 
form distribution and with a position uniformly distributed within 
±0.5 px in X and y of the SN location. The first approach involved 
repeating the original detection and aperture and PSF photometry 
routine on the pre-explosion images, with artificial stars inserted, 
and considering a detection to be any recovery of a star within 1 
pixel and 0.5 magnitudes of the input star's parameters. The sec- 
ond approach utilised isis to conduct image subtraction between 
the late-time images and pre-explosion images, in which artificial 
stars had been inserted in the latter. For this latter approach, we set 
the coordinates at which isis was to conduct aperture photometry 
and classified a detection to be any instance in which the recov- 
ered flux was 3 times that of the corresponding noise (including the 
systematic uncertainty; sec section 2.6). 

We consider the detection threshold to be the magnitude at 
which we recover 50% of input artificial stars, using a 3cr detection 
threshold with DAOphot. As noted by Maund (2013, in prep.), the 
completeness function can be considered in terms of the comple- 
mentary cimiulative Gaussian distribution. We therefore quote the 
corresponding width of the completeness function as an effective 
uncertainty on the derived detection threshold. 

3 OBSERVATIONAL RESULTS 
3.1 SN 1999ev 

SN 1999ev was discovered by T. Boles (Hurst et al. 1999) on 1999 
Nov 7.225 in the galaxy NGC 4274. Gamavich et al. (1999) sub- 



sequently classified the SN as being of Type 11, although no fur- 
ther sub-classification of the SN has been reported. Van Dyk et al. 
(2003a) attempted to identify the progenitor object in pre-explosion 
WFPC2 F555W images from 1995 Feb 1, although were not able 
to conclusively identify a single object as the progenitor In an in- 
dependent analysis, using a differential astrometric solution derived 
using post-explosion ACS WFC images containing the SN, Maund 
& Smartt (2005) were able to identify a star in the pre-explosion im- 
ages coincident with the SN position (with an uncertainty of 0.02"). 

Late-time ACS WFC F435VK, F555W and F814W images 
(with pixel scale 0.035" px ') of the site of SN 1999ev were ac- 
quired on 2010 Nov 14 (11 years post-discovery). A late-time im- 
age of the site of SN 1999ev is shown as Fig. 2. A geometric trans- 
formation was calculated between the post-explosion and late-time 
F555W images using 24 commons stars, with an uncertainty on 
the transformation of Ar = 0.016". A source is recovered in the 
late-time images at the transformed position of the SN as identi- 
fied in the post-explosion images by Maund & Smartt (2005), as 
shown in Fig. 3. The source is detected in all three filters: mf435w = 
25.63 ± 0.07, mF555w = 24.76 ± 0.06 and /Mf814w = 23.65 ± 0.05. 
The photometry of the SN in the post-explosion images was recal- 
culated, and the SN was measured to have mp^-jsw = 24.79 ± 0.10, 
rrifsssw = 24.19 ± 0.14 and mfsi4w = 23.49 ± 0.10. We note that 
the new measurement of the post-explosion F555 W photometry re- 
ported here is slightly fainter than measured previously, although 
the F435W and FS14W magnitudes are approximately similar to 
those of Maund & Smartt (2005). The 8.98 years between the post- 
explosion and late-time images reveals significant evolution in the 
light echo discovered by Maund & Smartt (2005), which has ex- 
panded to a radius of 0.48" from 0.25" (as shown on Figs. 3 and 
4). 

A transformation was calculated between the post-explosion 
and pre-explosion F5551V images using 24 stars (with a transfor- 
mation uncertainty of Ar = 0.039"). In the pre-explosion images, 
we identify the same source that Maund & Smartt (2005) identi- 
fied as the progenitor source with nifsssw = 24.66 ± 0.17. In ad- 
dition, we also find a nearby source with nifsssw = 25.08 ± 0.24 
located 0.2" (2 WF pixels) from the progenitor. We note that, for 
the period in which the pre-explosion observations were conducted 
the nearest logged warm pixel is 5 pixels away from the SN po- 
sition and not coincident with either the progenitor candidate or 
the nearby object. In the late-time image, however, we do not re- 
cover any source at the corresponding transformed position. The 
position of the pre-explosion source at the SN position was esti- 
mated using the three centring algorithms available to DAOphot 
(centroid, Gaussian and optimal filter) and the position determined 
using HSTphot PSF fitting. The positions are shown, with respect 
to the transformed position of the SN on the pre-explosion image, 
in Fig. 5. Although there is an apparent discrepancy in the posi- 
tions for the source and the transformed SN position, the discrep- 
ancy is not significant. It does, however, raise concerns about how 
positions are determined on subsampled images such as this pre- 
explosion WFPC2 WF2 F5551V image. The position determined 
using the optimal filter centring algorithm is noticeably different 
from the other three positions derived from the pre-explosion im- 
age and is offset in the direction of the nearby apparent neighbour- 
ing star. The standard deviation of the four measurements made on 
the pre-explosion image is 0.022". Given the apparent brightness 



* http://www-int.stsci.edu/ftp/instrument_newsAVFPC2/Wfpc2_hotpix/1995/ 
vary_9501 13_95021 1 J2.dat.Z 
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SN 1999evLATE 



Figure 2. Late-time colour image (composed of F435W, F5551V and 
FSWW images) of tlie area of NGC 4274 containing SN 1999ev. The image 
has dimensions 15" X 15", and is oriented such that North is up, and East 
is left. The position of SN 1999ev is located at the centre of the image. 




L ^ ^ 

Figure 4. Difference image between the post-explosion and late-time 
F5551V observations of the site of SN 1999ev. The evolution of the ex- 
panding light echo is evident; the white echo appears in the post-explosion 
image, while the outer dark echo occurs in the late-time frame. 

of the source, however, the astrometric uncertainty for the position 
derived using HSTphot alone may be as large as 0.04" or 0.4 WF 
pixels (Dolphin 2000b). 

3.2 SN 2003gd 

SN 2003gd was discovered by R. Evans on 2003 Jun 12.82, in 
the galaxy M74 (Evans & McNaught 2003). Kotak et al. (2003) 




Figure 5. Close-up of the pre-explosion site of SN 1999ev. In red (x) is the 
transformed position of the SN, derived from the post-explosion images, 
and the corresponding rm.s. eiTor ellipse. Also shown (in green) are the four 
separate measures of the position of the pre-explosion source derived using 
the centroid (□), Gaussian (A) and optimal filter (o) centring algorithms in 
DAOphot and the position determined using HSTphot (*). Also shown are 
the average (x) and the corresponding standard deviation ellipse of the four 
measurements. 

spectroscopically classified 2003gd as being a Type II SN, approx- 
imately 2 days post-explosion. Subsequent photometric and spec- 
troscopic observations of SN 2003gd, however, showed it to be a 
Type IIP SN discovered at the end of the plateau phase (Hendry 
et al. 2005). Smartt et al. (2003) made a preliminary identification 
of the progenitor in pre-explosion HS T WFPC2 F6Q6W and Gem- 
ini GMOS-N images. As such, SN 2003gd was the third SN, after 
SNe 1987 A and 1993 J, to have a progenitor identified in fortuitous 
pre-explosion images. Independent analyses by Smartt et al. (2004) 
and Van Dyk et al. (2003b) showed the candidate progenitor to be 
a RSG, corresponding to a star with initial mass of ~ 8 - 9Mq. The 
confirmation of this star as the progenitor was finally provided in 
2009, when the star was observed to no longer be present in late- 
time Gemini GMOS-N /' images (Maund & Smartt 2009); making 
it the first conclusively confirmed RSG progenitor for a Type IIP 
SN. 

Late-time ACS WFC observations of the site of SN 2003gd 
were acquired on 2010 Nov 14, and are presented on Figure 6. The 
SN position in the late-time images was determined with respect 
to the ACS HRC post-explosion images, with an uncertainty on 
the transformation of 0.010". The position of the SN on the pre- 
explosion HST and Gemini images, as presented by Smartt et al. 
(2004) and Maund & Smartt (2009), were recalculated to within 
0.028" and 0.025" respectively. The pre-explosion, post-explosion 
and late-time data F555W images are shown on Figure 7 and, for 
completeness, we also show the corresponding /' data presented by 
Maund & Smartt (2009) as Figure 8. 

In the late-time images we observe a source, termed Source 
A', that is clearly recovered in all late- time HST images at the trans- 
formed SN position, with magnitudes 25.93±0.04, 25.42±0.05 and 
24.90±0.04 in the F4351V, F5551V and F814W filters respectively. 
This source was detected in late-time Gemini GMOS-N g' and r' 
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Figure 3. HST imaging of the site of SN 1999ev. From left to right: Pre-explosion WFPC2 WF2 F555W image; Post-explosion ACS/WFC f 5551V image; 
and late-time ACSAVFC F555W image. 




images, but not recovered significantly in tiie corresponding im- 
age, acquired on 2008 Sep 06. Maund & Smartt (2009) measured 
g' = 25. 10+0.04 and r' = 24.49±0.05 (in Vega magnitudes) for the 
source at the SN position, and placed a detection limit of /' > 25.9. 
This source was also observed with HS T WFPC2 on 2007 Aug 1 1 
(for program GO - 11229; PI: M. Meixner). Photometry of these 
images using the HSTphot package yielded mf(,22w = 24.52 ± 0.08 
and infsi4w = 25.22 ± 0.26. We note that this photometry is ap- 
proximately 0.4 magnitudes brighter than the photometry of the 
same images reported by Otsuka et al. (2012). We confirmed our 
photometry using ISIS, determining the flux difference between the 
source in the 2007 WFPC2 F814W image and our late-time ACS 
FSMW image is consistent with the photometry conducted on the 
images directly. 

We recalculated the photometry of the source at the SN posi- 
tion in the pre-explosion WFPC2 F606W image, labeled Source A 
by Smartt et al. (2004), using HSTphot finding mfeoeiv = 25.06 ± 
0.06. Maund & Smartt (2009) derived the magnitude of the 
progenitor, using image subtraction techniques (see Figure 8), of 
23.85 ± 0.04 (with a possible 0.15 magnitude systematic uncer- 
tainty on underlying residual flux in the late-time Gemini image). 
Unlike the obviously red Source A observed in the pre-explosion 
HS T WFPC2 and Gemini GMOS images, it is apparent that Source 
A' in the late-time images is a blue-yellow object (see Figure 6). 
Even taking into account a colour correction between the late-time 
F814W and Gemini GMOS photometry (see Fig. 1), a signifi- 
cant increase in brightness is evident between the two observations 
separated by two years. 

3.3 SN 2004A 

SN 2004A was discovered by K. Itagaki (Nakano et al. 2004) on 
2004 Jan 9.4 in the galaxy NGC 6207. Kawakita et al. (2004) spec- 
troscopically classified the SN as a being a young Type II SN. 
Hendry et al. (2006) presented photometric and spectroscopic ob- 
servations of SN 2004A and showed it to be consistent with other 
normal Type IIP SNe, such as SN 1999em. Hendry et al. also pre- 
sented an analysis of the pre-explosion HS T WFPC2 observations 
of the site of SN 2004A from 2001 Jul 02, in conjunction with 
post-explosion ACS WFC observations of the SN acquired on 2004 
Sep 23. A source was barely recovered at 4.7cr at the SN position 
in the pre-explosion F8141V image. There was no corresponding 
source in the pre-explosion F6061V image, consistent with a star 



Figure 6. Late-time colom' image (composed of F435W, ^5551^ and 
FH14W images) of the area of M74 around the site of SN 2003gd. The 
image has dimensions 8" X 8", and is oriented such that North is up, and 
East is left. The position of SN 2003gd is located at the centre of the image. 



with ^6061^ - F814VK > 1.05. Hendry et al. concluded that if this 
was the progenitor star an RSG with initial mass 9^2^0' although 
given concerns about the significance of the detection of the source 
at the SN position, Hendry et al. placed a conservative limit on the 
initial mass of an undetected progenitor of < 12Mq. 

Late-time observations of the site of SN 2004A were acquired 
on 2010 Sep 09, approximately 6.7 years post-discovery. The 
late-time F8141V observation and the corresponding pre-explosion 
F814W observation are presented on Figure 9. The pre-explosion 
and late-time observations of the site of SN 2004A are presented 
in Figure 9. Using the post-explosion F8141V image, acquired on 
2004 Sep 23 with ACS WFC, the position of the SN on the pre- 
explosion and late-time frames was determined to within 0.021" 
and 0.028", respectively. The SN is not detected significantly in any 
of the late-time ACS WFC images. The 3cr detection limits at the 
SN position, in the late-time images, were evaluated with artificial 
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star tests to be mF4isw = 27.65, mfsssw = 27.5 and mpguw = 26.85 
mags. These limits are consistent with the expected depth for these 
images predicted by the ACS imaging exposure time calculator 

In the pre-explosion F814W image, HSTphot finds a source 
within 0.034" of the transformed SN position with mFsww = 
24.48 ±0.19 mags, detected with a signal-to-noise ratio of 5.6 (as 
shown on Figure 9). The positional uncertainty is slightly larger 
than the formal Icr uncertainty of the geometric transformation 
alone. The ability of HSTphot to determine the position of objects 
of such brightness, however, is limited, such that the expected un- 
certainty on the position of the object on the pre-explosion F8141V 
image is Ar > 0.060" (Dolphin 2000b). We note that this is the 
same object identified by Hendry et al. (2006) as the possible can- 
didate progenitor, although they measured the source to be ~ 0.2 
magnitudes brighter with a larger photometric uncertainty (we note, 
also, they did not correct their reported photometry to an infinite 
aperture, such that the discrepancy is a further 0.05 magnitudes 
larger). We also find the two sources A and B identified by Hendry 
et al. (2006) are hot pixels The pixel immediately adjacent to 
the pixel hosting the majority of the progenitor candidate's flux 
is a warm pixel; however this pixel has a low dark current (with 



low variabihty) " and was corrected by the OTFR pipeline. We 
determined the 50% completeness level for 3cr detections for the 
pre-explosion images using artificial star tests conducting the HST- 
phot. Artificial stars were placed in a 5px radius around the trans- 
formed SN position. We find the corresponding detection limits to 
be m/jftosw = 26.25 ± 0.40 and lUfgi^n- = 25.10 ± 0.45 mags. 

The late-time F814W image was subtracted from the pre- 
explosion image, and the difference images is presented on Fig. 
9. We find residuals at the location of the progenitor source and 
the two hotpixels. The disappearance of the pre-explosion progen- 
itor candidate in these late-time images confirms the authentic- 
ity of the object as the progenitor We derived photometry of the 
residual at the transformed SN position of m/rgMw = 24.57 ± 0.12 
mags, including a systematic uncertainty of ~ 10%. Using the re- 
sults from the artificial star tests, we derived a CTI correction of 
-0.216 ± 0.013, yielding a final magnitude for the progenitor of 
niFiuw = 24.36 ± 0.12 mags. This is similar to the photometry de- 
rived by Hendry et al. (2006), although for very different reasons 
and improved precision. 

A similar difference image was determined for the pre- 
explosion F606W image and the late-time F555VK image. No sig- 
nificant residual was found in the difference image as expected. 



' littp://etc.stsci.edu/etc/input/acs/imaging/ 

http://www.stsci.edu/hst/wfpc2/analysis/wfpc2Jiotpix.html 



" http://www-int.stsci.edu/instruments/wfpc2/Wfpc2_hotpix/2001/ 
vary_010617_010711_3.dat.Z 
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Figure 9. The site of SN 2004A in NGC 6207. Left) Pre-explosion WFPC2 F8UW image from 2001 Jul 02 with scale 0.1" px"'. The centred cross haks 
mark the transformed position of SN 2004A and the nearby sources A and B are known hot pixels on the WFPC2 WF2 chip. Centre) Late-time ACS/WFC 
F81AW image. Right) Difference image between the late and pre-explosion f 814W images. The observed residual soitfces correspond to the two hot pixels, 
A and B, and the progenitor. 



given the absence of a source at the transformed SN position in the 
pre-explosion image. Due to differences between the pre-explosion 
F606W and late-time F555W filter transmission fimctions, we did 
not use the difference image to derive detection limits for the pro- 
genitor. 

3.4 SN2005CS 

SN 2005cs was discovered by Kloehr et al. (2005) on 2005 Jun 
27.933 in the galaxy M51. Modjaz et al. (2005) spectroscopically 
classified SN 2005cs as being a young Type II SN. Richmond & 
Modjaz (2005) provisionally identified a blue supergiant in the field 
as a possible candidate for the progenitor, although later analysis (in 
conjunction with high resolution post-explosion HST ACS HRC 
images) by Maund et al. (2005) and Li et al. (2006a) found the 
progenitor star to be a RSG with initial mass Mzams ~ 8Mq. 

Late-time observations of the site of SN 2005cs were acquired 
on 2010 Jul 30 with the ACS/WFC, 5.1 years post-discovery. A 
comparison of the pre-explosion and late-time observations of the 
site of SN 2005cs is shown as Figs. 10 and 1 1 . In this case, the late- 
time observations exactly match the pre-explosion observations, us- 
ing the same filters and detectors. The pre-explosion and late-time 
observations were drizzled to a final common pixel scale of 0.035". 
Utilising post-explosion ACS HRC observations of SN 2005cs, the 
SN position was located on the pre-explosion and late-time images 
to within 0.007" and 0.004", respectively. 

Direct photometry of the source detected at the SN position in 
the pre-explosion F8I41V yielded mfgHw = 23.382 ± 0.048, which 
is 0.1 magnitude fainter than reported by Maund et al. (2005), but 
~ 0.3 magnitudes brighter than reported by Li et al. (2006b). In 
the late-time images, we do not recover a source at the transformed 
SN position. The detection limits in these filters were probed us- 
ing artificial star tests, yielding nip^^^w - 24.30 ± 0.4, nifsssw = 
24.95 ± 0.45 and mfSMw = 24.55 ± 0.15 mags. These limits are 
particularly high, relative to the expected depth for ACS/WFC im- 
ages of these durations, due to extended emission from the nearby 
cluster overlapping SN position. 

As noted by Maund et al. (2005) and Li et al. (2006b), this un- 
derlying emission can complicate the determination of the photom- 
etry of the progenitor from the pre-explosion imaging alone. This 
highlights the importance of using image subtraction techniques 



to accurately derive the progenitor photometry (by subtracting the 
background emission that is constant at both epochs). The differ- 
ence image between the pre-explosion and late-time FS14W ob- 
servations is presented on Fig. 1 1 . We measure the brightness of 
the progenitor to be mfmw = 23.62 ± 0.07 mags, which is fainter 
than the brightness determined from direct photometry of the pre- 
explosion source (see above). This magnitude is also significantly 
fainter than the photometry of Maund et al. (2005), and slightly 
brighter than the photometry of Li et al. (2006a) (who attempted to 
account for pre-explosion flux at the SN position due to the nearby 
cluster). We note that we find no significant source at the SN posi- 
tion in the corresponding F435W and F5551V difference images. 

Artificial star tests, in conjunction with image subtraction 
techniques, were used to derive alternative detection fimits (see 
Section 2.7) for the pre-explosion F435W and F555W images. In 
the absence of a corresponding late-time F65SN ACS/WFC im- 
age, the detection limit on the pre-explosion F65&N image could 
only be derived using direct recovery of artificial stars on the pre- 
explosion frame. The photometric completeness functions for the 
pre-explosion observations in which the progenitor was not de- 
tected is shown on Fig. 12 and presented in Table 3. The detection 
of a residual in isis difference images requires only a significant de- 
gree of residual flux at the SN position and is less dependent on 
the amount of background flux than the direct recovery of artificial 
stars. There are differences between the detection limits derived on 
the pre-explosion images here and the limits presented by Maund 
et al. (2005) and Li et al. (2006b). These studies used combinations 
of the analytical noise expression and artificial star tests, and treated 
the effect of flux from the nearby cluster differently. We note that 
our detection limits derived using image subtraction techniques are 
significantly deeper, highlighting the importance of late-time im- 
ages even in cases where detections of the progenitor are dubious 
or unavailable. 

We also find that there are a number of other sources that are 
clearly variable between the pre-explosion and late-time images in 
the vicinity of SN 2005cs (as shown on Fig. 1 1). Inspection of the 
pre-explosion and late-time images shows that the apparent resid- 
uals in the difference images are associated with stars which are 
clearly brighter or fainter in the late-time images compared with 
the pre-explosion images. Given the density of stars in this field, 
compared with the sites of the other SNe considered here, it is to 
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Figure 10. Colour images of the site of SN 2005cs using ACSAVFC before 
explosion and at late-times. The position of the SN is indicated by the cross- 
hairs. In the pre-explosion image a red source is clearly visible at the SN 
position, and is found to be absent in the late-time image. Each image has 
dimension 6" X 6" , and is oriented such that North is up. East is left. 

be expected that there would be other variable sources in the field 
around SN 2005cs. The other residuals in the difference images 
are, therefore, consistent with other real variables, but the progeni- 
tor object is the only star to be absent in one of the two sets of the 
images. 

3.5 SN 2006my 

SN 2006my was discovered on 2006 Nov 8.82UT by K. Itagaki 
(Nakano & Itagaki 2006) in the galaxy NGC 4651. Stanishev & 
Nielsen (2006) spectroscopically classified the SN as being a Type 
II SN similar to SN 1999em. Li et al. (2007), Leonard et al. (2008) 
and Crockett et al. (2011) analysed the pre-explosion WFPC2 im- 
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Figure 12. Detection completeness functions for pre-explosion ACS WFC 
observations of the site of SN 2005cs for the F435W, F555W and F65SN 
filters (in which the progenitor was not detected). The dashed curves in- 
dicate the completeness functions derived from attempted direct recovery 
of artificial stars on the original pre-explosion images, while solid curves 
are for detection limits derived in conjunction with image subtraction tech- 
niques. 



Table 3. Detection limits for pre-explosion ACSAVFC F435W, F555W and 
F658N images of the site of SN 2005cs measured using artificial star tests 
and direct recovery and image subtraction techniques. Previously derived 
detections limits from Maund et al. (2005) and Li et al. (2006b) are shown 
for comparison 





Recovery 


ISIS 


Maund et al. 


Li et al. 


F435W 


24.35 ± 0.20 


25.8 ±0.1 


24.8 




F555W 


24.85 ± 0.30 


264 ±0.1 


25 


25.6 


F658N 


24.15 ±0.30 









ages of the site of SN 2006my. All three studies commented on the 
significant offset between the transformed SN position and a nearby 
source recovered in the pre-explosion F8 14W image. Leonard et al. 
and Crockett et al. concluded that the F8I4H' source was unrelated 
to the SN, and that the progenitor was not detected in either the 
pre-explosion F5551V or F814VK images. 

The pre-, post-explosion and late-time F555W and F814W 
imaging of the site of SN 2006my is shown on Fig. 13. Following 
the analyses presented by Li et al. (2007), Leonard et al. (2008) and 
Crockett et al. (201 1), we analysed the pre-explosion WFPC2 WF2 
F555 W and F814W covering the position of SN 2006my. The SN 
position, derived from post-explosion WFPC2 PCI images, was de- 
termined on the pre-explosion images using 19 common stars with 
a resulting uncertainty of 0.024"; larger than achieved by Crockett 
et al. (2011). The transformed position is found to be in the prox- 
imity of a cluster of bright pixels and, as previously found by Li 
et al. (2007), Leonard et al. (2008) and Crockett et al. (2011), the 
transformed position is not consistent with the position of the near- 
est source found by HSTphot in the pre-explosion F814W image 
(see Fig. 13). The source detected by HSTphot is located 0.97 pix- 
els from the transformed SN position, an offset significantly larger 
than the transformation uncertainty. The transformed SN position 
is close to the position of a source in the pre-explosion F555W im- 
age but, as noted by the previous studies, based on the sharpness 
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Figure 11. ACS WFC F814W observations of the site of SN 2005cs. Left) Pre-explosion F814W image (0.44 years prior to discovery). Tlie cross hairs indicate 
the transformed position of the SN and the identified progenitor candidate. Centre) Late-time F814W image (5.09 years post-discovery) with the progenitor 
candidate absent. Right) Difference image, between the pre-explosion and late-time F814W observations, clearly showing a residual at the SN position. 



value derived by HSTphot, it is not consistent with a point source 
and is located 1.19 pixels from the source detected in the F814W 
image. The shift between the pre-explosion F555W and F814W 
image was found to be very small: Ax = -0.0039 and Ay = 0.0395. 

Using HSTphot, the pre-explosion F814W source was mea- 
sured to have brightness 24.24 ± 0.18, which is approximately 
0.2 magnitudes brighter than found in the previous studies (which 
is due to the use of different HSTphot settings used here). Sec- 
tions of the late-time F555iy and F8141V images were trans- 
formed to match the pre-explosion images; resampling the ACS 
WFC 0.05" pixels to 0. 1". These images were processed using isis, 
and the resulting difference images are shown in Fig. 13. isis sig- 
nificantly detects a residual, at a distance of only 0. 11 pixels from 
the transformed position of the SN. Using 5 reference stars, and 
the DOLPHOT photometry derived for these stars in the late-time 
FSMiy image, we derive a magnitude of 24.88 ± 0.13 mags for 
the residual in the difference image. Artificial star tests on the pre- 
explosion images were used to derive the CTI correction, for a star 
at the position of the observed residual on the pre-explosion frame, 
of 0.016 ± 0.003. The residual in the F8I4W difference image 
corresponds, therefore, to an object with raf8i4M' = 24.86 ± 0.13. 
Given the coincidence of residual with the transformed SN po- 
sition, we conclude that the progenitor was detected in the pre- 
explosion F814W image. The fainter F814iy magnitude derived 
using ISIS, compared to our own HSTphot photometry of the pre- 
explosion images and the previously reported values, and the appar- 
ent discrepancy in the position of the pre-explosion source and the 
SN, most likely reflects that the source in the pre-explosion F814W 
image is a blend of the progenitor with a source due East of the SN 
position (which also skews the apparent position of the progenitor 
source in that direction). In the analysis of the pre-explosion and 
late-time F555W images no residual was found in the difference 
image. The possible nature of the pre-explosion F555W source is 
revealed in the late-time images as a complicated, extended back- 
ground feature. 

We used artificial star tests to probe the detection limit of the 
pre-explosion F555VK and F814H' images within a 10 pixel ra- 
dius of the SN position; deriving 50% completeness limits at 3cr 
of mF555W' = 26.15 ± 0.65 and mF8i4w = 25.05 ± 0.75 mags. 



4 ANALYSIS 

In previous studies (e.g. Maund & Smartt 2005; Smartt et al. 2004; 
Hendry et al. 2006; Li et al. 2006a; Maund et al. 2005), the photo- 
metric properties of the progenitor and surrounding stars were de- 
rived through comparison with the ideal supergiant colour sequence 
presented by Drilling & Landolt (2000). As noted by Maund (2013, 
in prep.) there are significant deficiencies with this approach; such 
as the requirement for colour transformation equations to transform 
the observed photometry to the photometric system of Drilling & 
Landolt. More recent studies (e.g. Van Dyk et al. 2012; Maund et al. 
2011; Van Dyk et al. 2011; Fraser et al. 2012; Maund et al. 2013) 
have shown the benefit in fitting directly to SEDs constructed from 
synthetic spectra with known parameters using the same filters as 
the observations. 

In considering the observed photometry of objects identified 
at the position of the target SN in the pre-explosion images we 
utilised the BIX nested sampling and BASIE Markov Chain Monte 
Carlo SED fitting packages described by Maund (2013, in prep.). 
These two packages allow us to comprehensively probe the effects 
of different stellar parameters on the interpretation of the observed 
photometry on a densely sampled grid of model stellar photometry, 
in the native filter system of the observations. By design, both of 
these SED fitting packages can handle detections and upper limits 
simultaneously, although for limited data (i.e. the number of de- 
tections is less than the number of free parameters) only the BASIE 
code can be used to explore the allowed parameter space rather than 
locate a unique solution. Crucially, we can explore the degenera- 
cies between the parameters (such as temperature and reddening), 
and implicitly account for correlations between the temperature and 
bolometric luminosity (through the bolometric correction). 

Here we use two families of stellar SED models: the ATLAS9 
(Castelli & Kurucz 2004) and MARCS (Gustafsson et al. 2008) 
models. Synthetic photometry of these models was conducted using 
our own codes. 

As we expect the progenitors to be cool RSGs (< 4500^"), 
we interpret the observed photometry (and upper limits) with re- 
spect to the 5Mo spherical MARCS SEDs, which have been suc- 
cessfully compared with observations of RSGs in a number of pre- 
vious studies (e.g. see Levesque et al. 2005; Davies et al. 2013). 
We assume that RSGs are well described by models with surface 
gravity logg = 0.0, and fit for the effective temperature (F^//) and 
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Figure 13. Pre-explosion and late-time HST images of the site of SN 2006my in the F814W (top row) and F555W (bottom row) bands. The transformed 
position of the SN and the positions of the sources detected in the pre-explosion F555 W and F81AW (labeled "source") images are indicated by the circles. 



the reddening. We consider the effects of foreground and host red- 
dening to be due to Galactic-hke dust (parameterised by E(B - V); 
Cardelli et al. 1989). To constrain the effect of reddening due to 
dust expected to be found around RSGs, we consider the redden- 
ing laws for graphite or silicate dust (parameterised by the optical 
depth Tv), contained in spherical shells with ratios for the inner and 
outer radii of Ro,„/Rm = 2 or 10, following Kochanek et al. (2012). 
For each SN site we adopt models with metallicities appropriate for 
that site (see Table 1). 

To provide an additional handle on the interstellar reddening 
towards each progenitor, we also consider the reddening towards 
the stars immediately surrounding the SN position. Due to the dis- 
tances of the host galaxies, the surrounding stars, which are se- 
lected based on the condition that they are detected in all three of 
the late-time images for each SN, are expected to be luminous OB 
stars. For each of the stars, the parameters T^ff and E(B- V) are de- 
rived with respect to the ATLAS9 models. We selected those mod- 
els from the ATLAS9 grid that are consistent with supergiant sur- 
face gravities (Laidler et al. 2008)'^. Total reddenings are derived 
with respect to a CardelU et al. (1989) Ry = 3.1 reddening law, ap- 
propriate for reddening and extinction due to interstellar dust. The 
SED fits were conducted with the BIX Nested Sampling package 
(as all stars, by selection, had three colour photometry), such that 
the Bayesian evidence could be used to exclude those stars with 
colours clearly inconsistent with OB stars. Furthermore, to consider 



http://www.stsci.edu/hst/HST_overview/ 
documents/synphot/AppA_Catalogs4.html 



the SEDs of compact clusters we adopt the model spectra produced 
using the starburst99 code (Leitherer et al. 1999). 

For stellar progenitors, we derive masses using our various lu- 
minosity estimates (depending on the type of dust and star) follow- 
ing the technique of Smartt et al. (2009). Smartt et al. use predicted 
luminosities for the end-phases of STARS stellar evolution mod- 
els (Eldridge & Tout 2004) to derive initial masses for progenitors 
with luminosities constraints derived from the observations. For a 
given luminosity, the progenitor is considered to lie in the mass 
range bounded at one end by the most massive star to end core He 
burning at that luminosity, and at the other by the least massive star 
to proceed to model termination (the onset of core Ne burning) at 
that luminosity (see Smartt et al. 2009, and their Fig. 1). We use 
STARS models calculated at integer initial masses, with the appro- 
priate metallicities, and interpolate to determine the luminosities at 
the end of core He burning and the begiiming of core Ne burning. 
This scheme characterises possible RSG progenitors. We also note 
that, at lower masses, some stars will undergo second dredge-up, 
causing them to become Asymptotic Giant Branch (AGB) stars that 
are cooler but more luminous than the similar mass stars that die as 
RSGs. Based on the observed temperature range of RSGs, derived 
using MARCS spectra (Levesque et al. 2005), we use a tempera- 
ture threshold of 3400^r, above and below which we consider stars 
to be RSGs or AGB stars, respectively. 

In deriving posterior probability density functions (pdfs) for 
the initial masses for the progenitors, we also consider the effect 
of prior information from the IMF. We apply a weighting factor to 
the posterior pdfs oc M"^ '' for a Salpeter (1955) IMF, to follow the 
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weighting scheme applied by Smartt et al. (2009) for their analysis 
of the T^ypt IIP SN progenitor population. 



4.1 SN 1999ev 

The presence of a source in the late-time images precludes the use 
of image subtraction techniques to further analyse the nature of the 
pre-explosion source. We note that the pre-explosion F555W mag- 
nitude of the source at the SN position is of similar magnitude to 
the source in the late-time F555W image. The difference between 
the pre-explosion and late-time photometry is not significant with 
a p-value of 0.71 (using a simple z-test). This lends support to the 
hypothesis that the source observed in the pre-explosion images at 
the SN position has been recovered in the late-time images, and that 
the original identification of the progenitor presented by Maund & 
Smartt (2005) is, at least partially, incorrect. We suggest three pos- 
sible scenarios for the nature of the source at the SN position: 

(i) The source at the SN position in the pre-explosion and late- 
time images is a host cluster that contained the now absent progen- 
itor. 

(ii) The source observed in the pre-explosion and late-time im- 
age is an unrelated star that is coincident with the line-of-sight to 
the SN. While the late-time observations do suggest a large, young 
stellar population hidden by the large dust sheet, the determination 
of the likelihood of a chance alignment is non-trivial. Given the 
astrometric coincidence < 0.04", it is likely to be very low. 

(iii) The source observed in the late-time images is an unre- 
solved light echo, and the pre-explosion source has now disap- 
peared. Given the observation of evolving light echoes around the 
position of the SN, and the obvious amount of dust in the vicinity 
of the SN, the apparent late-time brightness may be due to a light 
echo from dust immediately behind the SN. We find this scenario 
unlikely, as it requires the progenitor and light echo to have coinci- 
dentally similar brightness. 

Given the nature of the late-time three-colour imaging it is not pos- 
sible to unambiguously distinguish between the different scenarios. 
The late-time images were used to examine the consequences of the 
progenitor residing in a host cluster that was observed in the pre- 
explosion and late-time images. The shape of the source in the late- 
time images was measured using the ishape package (Larsen 1999). 
In each filter band, ishape returned a significantly better fit with a 
Moffat function over a delta function (;^^^(Mofl'at)/;('^ (delta) < 0.95) 
and an effective radius R^ff > O.lx the Full Width at Half Maxi- 
mum (Larsen 1999). The effective radius of the source was mea- 
sured to be 0.99:15^5, 1.82:5;^6 and 1.18!° pixels in the F435W, 
F555W and F814W, respectively (with a pixel scale corresponding 
to 2.7pcpx"' at the distance of NGC 4274). The large error bars 
are symptomatic of the complexity of the region hosting the SN, 
in particular with the proximity of light echoes and the apparent 
faintness of the source. The ishape analysis was also conducted on 
six nearby objects that were all found to be consistent with point- 
like, stellar sources - suggesting that ishape does have the capa- 
bility, under the conditions of the late-time images, to differenti- 
ate extended sources from point-like sources. To further explore 
the implications of a host cluster for the progenitor, the late-time 
ACS photometry was compared with Starburst99 models (Lei- 
therer et al. 1999); and the results of this fit is shown as Figure 
14. Given the three colour photometry, there are two allowed so- 
lutions: a moderately reddened older solution (40 - IOOM3T) im- 
plying MzAMs < 9Mo; and a heavily reddened younger solution 



(< lOMyr) implying Mzams > 20Mq'^. In addition, given the cri- 
terion presented by Bastian et al. (2005), that point-like objects 
with Mv < -8.6 are more likely to be clusters than individual 
bright stars, requires E(B - V) > 0.76 (assuming an Ry = 3.1 
Galactic reddening law). In terms of shape, absolute brightness and 
colours, the late-time source is consistent with a cluster. Using the 
photometry of nearby point-like sources, a weighted-average red- 
dening of £(B - V) = 0.95 ± 0.32 was measured. The large error 
bar is consistent with both the poor photometric errors for each 
of the six nearby sources and the large scatter in reddenings in 
the sample. The reddenings towards these objects are significantly 
greater than expected for just pure foreground Galactic reddening 
(E(B - V) = 0.2) towards NGC 4274. This may reflect a com- 
plex dust distribution where some of the sources are in front of the 
dust sheet and others may be embedded. This suggests that there 
might be significant reddening towards the source at the SN posi- 
tion, however this is not conclusive. 

If the late-time source is, in fact, a light echo, then the late- 
time images do not provide any further insight into the properties of 
the source in the pre-explosion images. Our inability to the confirm 
the disappearance of the progenitor also means that the late-time 
images cannot be used to rule out the possibility that the object in 
the pre-explosion and late-time images is an unrelated object in the 
line-of-sight. High-resolution near-infrared observations, with the 
HST, could be used to probe the nature of the stellar population 
behind the dust sheet (to examine the density of objects along the 
line-of-sight) as well as provide further constraints on the nature of 
the source in late-time images and nature of the obscuring dust. 
The ambiguity of the nature of the object at the SN position 
in the pre-explosion and late-time images means that, although 
SN 1999ev may have had an identified progenitor of some kind, 
the previously derived progenitor properties are unreliable; even in 
the interpretation that the source is a cluster, the reliance on three- 
colour photometry leads to degeneracies in the reddening-age solu- 
tions that prohibit a precise initial mass estimate for the progenitor. 

4.2 SN2003gd 

The photometry of 30 stars within 4" (~ 200pc) of the position of 
SN 2003gd was used to derive an average reddening towards the 
SN site of E(B - V) = 0.14 ± 0.04 (see Figure 15). The amount 
of reddening is consistent with the reddening previously estimated 
from three colour photometry of the surrounding stars using early 
post-explosion ACS HRC images (Smartt et al. 2004) and from the 
colour evolution of the SN itself (Hendry et al. 2005). 

Given the detection of the progenitor in pre-explosion obser- 
vations in two filters, we consider the roles of three different types 
of reddening: 1) unconstrained reddening, with an interstellar red- 
dening law: 2) an interstellar reddening component consistent with 
the observed reddening to the surrounding stars and an uncon- 
sttained degree of reddening arising from Graphite dust around the 
progenitor; and 3) the same as 2, except with Silicate dust. The 
corresponding regions of the parameter space and the Hertzsprung- 
Russell (HR) diagram, allowed by the pre-explosion observations 
in conjunction with half-solar metalhcity MARCS SEDs, are pre- 
sented on Fig. 16. Given the observed colour of the progenitor, re- 
gardless of the amount of reddening, we find T^ff > 3500^" and the 

We note that the starburst99 code uses the Padova stellar evolution (Gi- 
rardi et al. 2002) models, and the masses we report for the derived ages are 
from these models. 
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Figure 15. Reddenings for stars surrounding the sites of SNe 2003gd, 2004A, 2005cs and 2006my (as a function of distance from the SN positions). 





E(B-V) 

Figure 14. The posterior probability distribution for SED fits of the late- 
time photometry of the source at the position of SN 1999ev to cluster SED 
models produced using Starburst99. The contours contain the 68% and 
95% of the probability, while the most likely solution is indicated with the 
cross. The vertical dashed line corresponds to E(B -V) = 0.76; for E{B - 
V) > 0.76, the absolute magnitude of the source is Mps^syf < -8.6. The 
horizontal lines indicate the lifetimes of stars with a given initial mass as 
predicted by the Padova stellar evolution code (Girardi et al. 2002). 



E{B - V) = 0.06 ± 0.03, derived using three colour photometry 
from the post-explosion ACS IWFC images. 

Due to the single F814VK detection of the progenitor, and 
the upper ^606^ limit, an SED fit with > 1 parameter is under- 
constrained by the observations. Using particular prior assump- 
tions, such as the reddening derived from surrounding stars, it is 
possible to probe the allowed likely parameter space. We con- 
ducted a two-parameter fit to MARCS SEDs, with half-solar metal- 
licity, for temperature and a single reddening component (follow- 
ing a Cardelli et al. reddening law) with the reddening derived 
from surrounding stars as a prior. The resulting constraint in the 
temperature-reddening plane is shown on Figure 17. The low red- 
dening and the strict F606W upper limit constrains the temperature 
of the progenitor to be lower than 3700K; although this is sensitive 
to the reddening prior, such that higher reddenings may permit hot- 
ter progenitors. The corresponding region on the HR diagram is 
also shown on Figure 17. As the reddening is effectively fixed, the 
slope of the contours reflects the increasing bolometric correction 
for cooler stars. If the reddening inferred from the surrounding stars 
is, instead, a lower limit on the reddening towards the progenitor 
(due to an additional component of reddening due to circumsteller 
dust), then we can only derive a lower limit on the luminosity of 
the progenitor. The contours also include those stars in the initial 
mass range 5-7Mo that are expected, at this metallicity, to undergo 
second dredge up and become AGB stars. 



radius of the progenitor constrained to be 200 < R < 400./?q. The 
flatter nature of the graphite and silicate reddening laws, compared 
to the Cardelli et al. reddening law, leads to tighter luminosity con- 
straints and a weaker dependence on the temperature/colour of the 
progenitor. The masses inferred for the progenitor are relatively in- 
sensitive to the choice of reddening law, with M,„„ ~ 7 - 8 ± I.QMq. 

Given the reddening derived towards SN 2003gd, we find that 
the SED of the object recovered at the SN position in the late-time 
images (A') is consistent with a black body with temperature T ~ 
9 OOOK, which is similar to the expected temperature for SNe at 
such late times (see e.g Otsuka et al. 2012). 



4.3 SN2004A 

The three-colour photometry of 12 stars within 5.25" (~ 500pc; sec 
Fig. 15) of the position of SN2004A was used to derived a weighted 
average reddening towards the SN site of E(B - V) = 0.16 ± 0.06. 
This is larger than the estimate made by Hendry et al. (2006) of 



4.4 SN 2005CS 

The photometry, from the pre-explosion images, of 20 stars within 
2" (~ 120pc) of the position of SN 2005cs were used to determine 
a weighted average reddening towards the SN of £(S - F) = 0.22 ± 
0.05 (see Fig. 15). 

In order to constrain the properties of the progenitors we 
utilised the F8141V magnitude presented here, measured using im- 
age subtraction, and the revised detection thresholds for the pre- 
explosion ACS F435 W, F555 W and F65SN images. Furthermore, 
we adopted the infrared detection limits for the progenitor reported 
for pre-explosion Gemini NIRI JHK images (Maund et al. 2005) 
and NICMOS FllOW, F160W F222M images (Li et al. 2006a). 
As the NICMOS images are deeper, they form the principal con- 
straint on the SED of the progenitor in the IR, however we include 
the Gemini NIRI limits in our calculation for completeness. We 
explored the same parameter space, for the same combinations of 
reddening components, as for the progenitor of SN 2003gd (see 
Section 4.2) and the results are presented on Fig. 18. 

Similarly to the progenitor of SN 2003gd, the graphite and 
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Figure 16. The parameters of the progenitor of SN 2003gd for different reddening components. In each panel the contours contain 68% and 95% of the 
total probability. Top Row), Left) The temperature and reddening of the progenitor assuming Galactic-like dust; Centre) the progenitor's location on the HR 
diagram (also shown are the locations of stellar evolution models of given initial mass for the end of core He burning (A), the onset of Ne burning (■) and the 
endpoints for those models that undergo second dredge-up (•). Dotted grey lines indicate lines of constant progenitor radius.; and Right) the inferred initial 
mass probability density function with no weighting (solid Une) and weighting according to the initial mass function (dotted line). Middle Row) The same as 
the top row, but for a Cardelli et al. (1989) reddening component of E(B - V) = 0. 14 ± 0.04, derived from the surrounding stars and a unconstrained reddening 
component for graphite dust around the progenitor (solid contours are for Rout/^in = 2 and dotted contours are for R„utlRin = 10). The mass probability 
density function are as for the top row, but only shown for the Routl^in = 2 solution. Bottom Row) The same as the middle row, but for silicate dust. 



silicate reddening laws yield flatter contours of the HR diagram, 
making the dependence of the luminosity on the temperature less 
extreme. For each reddening tjqje, there are two islands of preferred 
solutions: a cool, low reddened solution and a hotter, reddened so- 
lution; reflecting the severe constraints from non-detections in the 
infrared and in blue, respectively. The bimodal probability distribu- 
tion in Tgff and E{B — V) leads to a skewed mass probability den- 
sity function for a Cardelli et al. reddening law extending to higher 
masses, although the peak of the distribution occurs at ~ lOMg. 
For the graphite and silicate reddening laws, the unweighted mass 
probability density function is more symmetric and peaks around 
~ UMq. The strict infrared Umits exclude the possibility of the 



progenitor being a massive AGB stars for all the reddening types 
(Eldridge et al. 2007). 



4.5 SN 2006my 

Previously Li et al. (2007) suggested that, apart from a Galactic 
reddening component, there was no evidence for a host redden- 
ing component in spectra of SN 2006my; this value was similarly 
used by progenitor studies conducted by Crockett et al. (201 1) and 
Leonard et al. (2008). We utilised the HSTphot photometry of the 
post-explosion WFPC2 images of SN 2006my to study the redden- 
ing associated with the surrounding stellar population. We selected 
46 good stars ix^ < 1.5, |sharp| < 0.3) with complete three-colour 
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Figure 17. The same as Fig. 16 but for pre-explosion observations of the progenitor of SN 2004A, assuming only a Cardelli et al. reddening law with 
£■(8 - y) = 0.06 ± 0.03 derived from surrounding stars. 
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Figure 18. The same as Fig. 16, but for the pre-explosion observations of the progenitor of SN 2005cs. 
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photometry within 5" (~ 540pc) of the SN. This photometry was 
compared with solar metallicity supergiant SEDs, to derive an av- 
erage reddening towards SN 2006my of E{B- V) = 0.49+0.25 (see 
Fig. 15). The large scatter reflects the poor quality of the WFPC2 
photometry, as well as likely differences in the reddening between 
the individual stars. The relative colour constraint provided by the 
pre-explosion F555W limit and the F814W detection of the pro- 
genitor is insufiicient to limit the temperature of the progenitor for 
temperatures below < 4500^r. The allowed regions of the HR di- 
agram, for the two reddening estimates, is shown on Figure 4.5. 
As for SN 2004A (see Section 4.3), the shape of the contours is 
dictated by the larger bolometric correction at cooler temperatures. 
The lack of an infrared detection or detection limit for the progen- 
itor means that both the RSG and AGB solutions are allowed. The 
higher reddening inferred from the photometry of the surrounding 
stars implies higher luminosities for the progenitor, than for just 
Galactic foreground reddening, and leads to a higher initial mass 
(13.4 ±2.8 vs. 9.8+ 1.7Mo). 

Given the apparent offset between the transformed SN po- 
sition and the source recovered in the pre-explosion F814Ty im- 
age, we explored the possible causes for the apparent discrepancy. 
We note that all three previous studies, and our own results, agree 
that there is a significant offset between the transformed SN po- 
sition and the pre-explosion ^814^ source. As noted by Leonard 
et al. (2008), in conjunction with Dolphin (2000b), the positional 
uncertainty for an isolated source with the brightness of the pre- 
explosion FS14W source is ~ 0.4 pixels. We conducted Monte 
Carlo simulations, where the pre-explosion image was "resampled" 
under the assumption that the observed counts could be modelled as 
source and background flux (following a Poisson distribution) with 
a Gaussian readnoise contribution. Aperture photometry of these 
simulated images was conducted using DAOphot, with the ofilter 
centring algorithm, to recover the position of the source (within a 
5 pixel centring box). The result of these simulations is shown on 
Fig. 20 and the average position isx = 410.5 ± 0.5, y = 159.2 ± 0.6; 
offset from the SN position by 0.4 pixels. From the outcome of the 
Monte Carlo simulations we note three effects: 

(i) The positional uncertainty is larger for faint sources and will 
be dependent on pixel noise statistics in the main pixel containing 
source flux and the surrounding pixels. 

(ii) Given the subsampled nature of the PSF in WFPC2 images, 
the uncertainty on the position of a source will also be dependent 
on the brightness of the immediately surrounding pixels; and, as 
evident from Fig. 20, the measured position may be skewed by the 
proximity of the true source position to a pixel edge. 

(iii) The position derived is sensitive to the choice of centring 
algorithm utilised. Trials using the "centroid" centring algorithm 
in DAOphot showed that positions were preferentially located at 
integer and half-integer pixel coordinates. 

We conclude that, in addition to uncertainties associated with de- 
termining positions on WFPC2 images for isolated sources, it is 
also important to consider the eflfects of nearby sources (within a 
few pixels) that may skew/bias the centring algorithm away from 
the true source position. From the Monte Carlo simulations, given 
the environment and pixel noise at the SN position, we estimate 
the imcertainty on the position of the pre-explosion FS14W source 
may be ~ 0.078". With such large uncertainties, the apparent off'set 
between the SN and the pre-explosion source would be only ~ Icr 
given the offsets calculated by Crockett et al. (2011) and Leonard 
et al. (2008); however, we caution that our lower quality geometric 
transformation results in a 1.2cr off'set. 



5 DISCUSSION 

A summary of the masses derived for the progenitors considered 
here is presented on Table 4. The biggest dififerences between our 
findings and those of Smartt et al. (2009) are for the progenitors of 
SNe 1999ev and 2006my. Our late-imaging shows that the nature 
of the source found at the SN position in pre-explosion observa- 
tions of the site of SN 1999ev is, at best, uncertain; whilst, at worst, 
a misidentification of a host cluster or unrelated coincident star. As 
the nature of this source is not clarified in the late-time imaging, the 
progenitor of SN 1999ev should no longer be considered in progen- 
itor population statistics. In the analysis presented by Smartt et al. 
(2009), the progenitor of SN 1999ev had the distinction of having 
the highest mass inferred for a detected Type IIP SN progenitor, 
and so set the maximum mass limit for stars to explode as RSGs 
and produce Type IIP SNe. With the removal of the SN 1999ev 
progenitor from the population statistics, the maximum mass limit 
will actuaUy drop and make the "Red Supergiant Problem" even 
more severe. The late-time imaging of the site of SN 2006my has 
shown that the progenitor was detected in the pre-explosion obser- 
vations, and the upper mass limit quoted by Smartt et al. should 
now be quoted as a detected progenitor with a corresponding mass 
estimate. 

The mass estimates derived for the confirmed progenitors are 
generally higher (by ~ IMg) than those presented by Smartt et al.; 
in part, due to the slightly larger foreground and host reddenings we 
inferred towards the progenitors from the colours of the surround- 
ing stars. In addition, for SN 2003gd and 2005cs we considered 
additional reddening components to the reddening from interstel- 
lar dust, in keeping with the expectation that there is dust local to 
progenitor that is not probed by the surrounding stars or the obser- 
vations of the SNe. We also note that our uncertainties are smaller 
than those quoted by Smartt et al.: we considered the effects of 
uncertainties of the luminosity convolved with the flat probabil- 
ity distribution of the star having a mass in the range bounded by 
the maximiun mass star to end core He-burning at that luminosity 
and the minimum mass star to begin Ne burning at that luminos- 
ity. The initial mass probability density functions are, apart from 
SN 2005cs, approximately symmetric and almost follow a normal 
distribution. We believe this is a fairer presentation of the initial 
masses for the progenitors and their uncertainties. The application 
of a weighting to the initial mass pdf, according to the IMF, has 
a small effect in shifting the pdf to slightly lower masses. Unlike 
Smartt et al. (2009), who used this weighting scheme to "truncate" 
their large uncertainties, we find that the efifect of weighting on the 
relative width of the initial mass pdfs is minor. 

The sUght increase in the inferred progenitor mass does not 
help rectify the apparent discrepancy between these "evolutionary" 
masses and the progenitor masses derived from hydrodynamical 
models. Utrobin & Chugai (2008) found an initial mass for the pro- 
genitor of SN 2005cs of 18.2 ± IMq, which is at odds with the 
new masses derived here, regardless of the choice of reddening. As 
Utrobin & Chugai (2008) suggest, additional dust in a circumstel- 
lar shell could lower the apparent luminosity of the progenitor and 
decrease the mass inferred from pre-explosion observations. The 
nature of the pre-explosion observations of SN 2005cs, in particu- 
lar the strict near-infrared upper limits, severely limits the amount 
of reddening that the progenitor might undergo. The constraints on 
the radius for the progenitor of SN 2005cs are also below the ra- 
dius inferred by Utrobin & Chugai (2008) of 600 ± UOR^, but not 
significantly discrepant. 

It is clear from the analysis presented in Section 4, that the 
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Figure 19. The parameters of the progenitor of SN 2006my. Left) The progenitor of SN 2006my on the HR diagram, assuming E{B — V) = 0.49 ± 0.25 
(solid contours) and E(B ~ V) = 0.027 (dotted contours). Right) The initial mass probability density functions for the progenitor of SN 2006my, given 
E(B -V) = 0.49 ± 0.25 (heavy curves) and E(B -V) = 0.027 (grey curves). 
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Figure 20. The results of position determinations for the pre-explosion source, using DAOphot and the ofilter centring algorithm, for 1000 Monte Carlo 
simulations of the pre-explosion f 8 14W image of the site of SN 2006my. Left) The positions of the recovered sources in the coordinate frame of the original 
WFPC2 image (in the DAOphot coordinate system). The position of the source recovered by HSTphot is indicated by the +, while the mean position derived 
form the Monte Carlo simulations is indicated by the * ; Centre and Right) The position distributions in x and y coordinates for the recovered sources; the grey 
line indicates a Gaussian fit to the distributions. 



available pre-explosion observations of a given SN progenitor dic- 
tates the degree of analysis that may be conducted. Given the pres- 
ence of constraining optical and infrared upper limits, the possible 
effect of circumstellar reddening on the progenitor of SN 2005cs 
was evaluated despite having a detection in only one band. This 
demonstrates the importance for having good pre-explosion ob- 
servations in the near-infrared for studying the cool progenitors 
of Type IIP SNe, even if the progenitor is not detected at those 
wavelengths. The two detections of the progenitor of SN 2003gd, 
at different wavelengths, enabled similar constraints for reddening 
due to circumstellar dust. Conversely, the analyses of SN 2004A 
and 2006my were limited by them having only a single detection 
and only loose constraint on the progenitor colour at a bluer wave- 
length. 

Due to the lack of constraints on reddening due to circum- 
stellar dust, the final derived masses for the progenitors of SNe 
2004A and 2006my may represent, in actuality, lower mass lim- 
its. For this study, we have assumed the circumstellar dust follows 
the reddenings laws proposed by Kochanek et al. (2012), for the 
case of the progenitor of SN 2012aw. Conversely, Van Dyk et al. 



(2012) suggested a steeper reddening law, still following Cardelli 
et al., but with Ry ~ 4.35; leading to an overall higher extinction. 
The Kochanek et al. (2012) reddening law was determined using 
models of specific dust compositions, of either graphite or silicate 
dust, expected to be found around RSGs, whereas Van Dyk et al. 
(2012) estimated the change in reddening law based on observa- 
tions of Galactic RSGs. For SN 2012aw (Fraser et al. 2012; Van 
Dyk et al. 2012), the progenitor was detected in four bands, and 
the degeneracies between reddening, reddening law and tempera- 
ture could not be broken; suggesting the full determination of the 
parameters, independent of the assumptions of reddening laws, re- 
quires detections at > 4 wavelengths, such as for the progenitor of 
SN 2008bk (Mattila et al. 2008; Van Dyk et al. 2012). 

Both Fraser et al. (2012) and Van Dyk et al. (2012) noted 
the large decrease in reddening determined for the pre-explosion 
source and for the subsequent SN; suggesting such large redden- 
ings may affect all progenitors but not be apparent post-explosion. 
For their sample of RSGs, Davies et al. (2013) measure extinctions 
arising from circumstellar dust in the range Ay = 0.0 - 1.0. A fur- 
ther issue, that we have not explored, was suggested by Walmswell 
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& Eldridge (2012, and references therein) that the amount of dust 
in the circumstellar medium is related to the mass loss rate of the 
RSG and, ultimately, its bolometric luminosity. 

A corresponding issue to the reddening problem is the tem- 
perature. Wc note that, for the progenitors with constraining pre- 
explosion observations, we find the allowed temperature range to 
be generally hotter than the predicted endpoints for the stellar evo- 
lution models, but are consistent with the recent reappraisal of RSG 
temperatures by Davies et al. (2013). The lower limit of the temper- 
ature scales for the progenitors of SNe 2003gd and 2005cs suggests 
that they have a spectral type no later than MO, which corresponds 
to the predicted positions for stars that have just finished core He- 
burning; the pre-explosion observations of SNe 2003gd and 2005cs 
suggest that the progenitors were not massive AGB stars (Eldridge 
et al. 2007; Siess 2007). It is only in the poorly constrained cases, 
for the progenitors of SNe 2004A and 2006my, that we cannot ex- 
clude cooler temperatures that might be associated with massive 
AGB stars. With limited observations in the optical (in particular 
the B and V bands), it is difficult to place limits on the maximum 
temperature of the progenitor, as hotter temperatures can always be 
acconmiodated with additional reddening. In the case of the pro- 
genitor of SN 2005cs, the requirement that the progenitor ended its 
life as an RSG has serious implications for the interpretation of the 
subsequent SN as a low luminosity "Electron-Capture" SN (Janka 
2012). 

Our late-time imaging campaign has shown that, for the case 
of SN 2003gd, the possibilities of recovering precise photometry 
of the progenitor through template subtraction may be undermined 
by rebrightening of the SN at late-times. In the case of SN 2003gd, 
the previous analysis of Maund & Smartt (2009) was fortuitous in 
that it managed to observe the SN before it rebrightened with a 
very strict brightness limit of the SN in the late-time (' image. Such 
rebrightening is not without precedent; Kotak et al. (2009) observed 
the optical lightcurve of SN 2004et to rebrighten (by ~ 1 mag in V) 
at optical and infrared wavelengths. 

As noted in Section 4.5, the apparent discrepancy between the 
position of the progenitor of SN 2006my on the pre-explosion im- 
ages and the transformed SN position does raise questions about 
how positional uncertainties are handled. Maund & Smartt (2005) 
determined the positional precision using the standard deviation of 
the four different centring methods available to DAOphot (centroid, 
ofilter, gauss and psf). It must be noted, however, that these centring 
techniques are not providing independent estimates of an object's 
location. In the future, it may be preferable to choose a single cen- 
tering technique, and consider the role of Poisson noise (both object 
and background) and read out noise in each pixel on the determina- 
tion of an object's position. Furthermore, tests using the centroiding 
routines in both the DAOphot and SAO image DS9 packages has 
shown that, in the case of isolated objects in subsampled images, 
centroiding can be a relatively blunt tool (providing default posi- 
tions located in the centre of pixels). In considering ffux deficits 
using image subtraction techniques, rather than appealing to as- 
trometric coincidence, we have shown that confirmation of a star 
as being the actual progenitor requires observing it to have disap- 
peared. 



6 CONCLUSIONS 

We have presented late-time imaging of the sites of five Type IIP 
SNe with pre-explosion HST images, in which progenitor candi- 
dates were detected. In three of the cases (2003gd, 2004A and 



2005cs), our previous identifications have been confirmed and we 
find initial masses for these stars in the range 6 - 14Mq. The pre- 
explosion observations of SNe 2003gd and 2005cs are sufficient to 
place constraints on the progenitor mass that are relatively insensi- 
tive to the amount and type of dust around these progenitors. Given 
the similarities in brightness between the pre-explosion and late- 
time sources detected at the position of SN 1999ev, we conclude 
the progenitor identification for this SN is unsafe and suggest the 
pre-explosion source may be a reddened host cluster; although the 
three-colour late-time imaging is insufficient to place a tight con- 
straint on the age or reddening of a such cluster. The analysis of the 
pre-explosion and late-time observations of the site of SN 2006my 
have revealed that the source previously thought to be significantly 
offset from the SN position has disappeared. The astrometric coin- 
cidence of the residual in the difference image with the transformed 
SN position suggests it was the progenitor object. 

Far from providing just simple confirmation of a progenitor's 
identity (through its disappearance), our analysis shows late-time 
imaging is crucial for conducting a deeper and more precise anal- 
ysis of the properties of a progenitor than is afforded by fortuitous 
pre-explosion observations alone. The power of the application of 
late-time imaging, for studying progenitors, is demonstrated by the 
significantly deeper detection limits that may be achieved by us- 
ing artificial star tests in conjunction with image subtraction tech- 
niques. 
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